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Equations for a perfect fluid can be obtained by means of the variational principle both 
in the Lagrangian description and in the Eulerian one. It is known that we need additional 
fields somehow to describe a rotational isentropic flow in the latter description. We give 
a simple explanation for these fields; they are introduced to fix both ends of a pathline in 
■^T the variational calculus. This restriction is imposed in the former description, and should be 

imposed in the latter description. It is also shown that we can derive a canonical Hamiltonian 
formulation for a perfect fluid by regarding the velocity field as the input in the framework 
of control theory. 



£h ' §1- Introduction 

>> 

The Euler equation, together with mass and entropy conservations, describes 
dynamics of the perfect fluid, and can be derived from the variational principle in 
the Lagrangian description. 1 - Let us write p for mass per unit volume, s for entropy 
per unit mass, and v for velocity field. Because of the local equilibrium, the internal- 
energy density per unit mass, e, is a function of p and s. Apart from constraints 

£>y coming from the conservation laws, the Lagrangian density is given by 

CU .1.2 



2 ^-e(p, S )| . (1-1) 

J> ■ The action, which is minimized to yield the Euler equation, is given by the integral 

of the Lagrangian density over the space and time considered. In the Lagrangian 
description, the position should be regarded as a variable of the Lagrangian density, 
and thus the velocity field v is given by the time derivative of the position of a fluid 

^. . particle. 

In the Eulerian description, not the position but the velocity is a variable of 
the Lagrangian density. Minimizing the action, we can obtain the Euler equation. 
However, unlike in the Lagrangian description, the resultant velocity field cannot 
be rotational on the isentropic condition. To remove this flaw, BatemaEp* added 

/\ ■ some scalar fields, sometimes called Clebsch potentials p*™ to the Lagrangian den- 

sity. Later LirP 1 used more additional scalar fields, which Selinger and WhithairP 1 
considered to be redundant. Schutz-^ also used Bateman's scalar fields in considering 
the general-relativistic gravitation field, and still complained of too many additional 
scalar fields. Irrespective of these controversies, the variational principles^ 1 '^'^ in 
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the Eulerian description have been used in some studies on magnetohydrodynam- 
ics, multivalued plasmas, elasticity and a hydrodynamic description of relativistic 
starsP*'® ™ Recently, Kambe^ 1 has explained Clebsch potentials from the point 
of view of symmetry of the gauge theory. Yoshida^'^ claimed that the number 
of required additional fields should be more than Bateman's and fewer than Lin's. 
Some authors^ 4s«>l§)il!2H15) considered that Clebsch potentials are related with the 
Lagrangian coordinates. 

We show below that Clebsch potentials in the variational calculus of the Eu- 
ler description are introduced to fix both ends of a pathline; the same restriction 
is imposed in that of the Lagrangian description. In relation to the variational 
principles, various canonical and noncanonical Hamiltonian formulations have been 
proposed.^"^ We derive a canonical Hamiltonian formulation by means of an 
optimal control theory known as Pontryagin's minimum principle. 

We state the problem in more detail and show our notation in |2j Our expla- 
nation for Clebsch potential is given in ^3) We give a brief review of the control 
theory in Q4.1\ and derive a canonical formulation by means of the control theory in 
£ 34.21 The last section is devoted to discussion. We compare our study with previous 
works on Clebsch potentials and Hamiltonian formulations in ^5j Some details are 
relegated to Appendices. Although we limit the following discussion in the text to 
nonrelativistic perfect fluid, we show in Appendix ICl that our discussion is also valid 
for a relativistic perfect fluid. 

§2. Statement of the problem 

We write V for the spacial region of the container filled with a perfect fluid, and 
dV for the surface of the container. We fix the container and consider the dynamics 
of the fluid from the initial time £; n it to the final time t^ a . Let r denote time in 
the Lagrangian description and t denote time in the Eulerian one although they 
are equivalent in the nonrelativistic theory. The partial derivative with respect to 
r (d r ), and t (dt), imply the Lagrangian and Eulerian time derivatives, respectively. 
See Appendix [A] for their relation. We label a fluid particle with its initial position 
a, and write X(a,r) for its position at time r. Thus, a = X(a, ti n it) gives the 
Lagrangian coordinates. The volume element in the Lagrangian coordinates can be 
given by the determinant of the Jacobian matrix, 

t(„ ^ _ d(Xi,X 2 ,X 3 ) 
d(a 1 ,a 2 ,a 3 ) 

where Xi and a, are, respectively, the components of X and a. By definition, we have 
J(a, iinit) = 1- We assume that a fluid particle never shrinks to a point, i.e., J(a,r) 
has no singular points, in the space and time considered. Thus, we can define the 
inverse of x = X(a,r), for which we write a(x,t). Let us write T for temperature 
and p for pressure, and the first law of thermodynamics, de = —pdp' 1 + Tds, yields 



where the subscripts s and p indicate variables fixed in the respective partial differ- 
entiations. Below, Roman indices run from 1 to 3 except in 2] and Appendix[Bl and 
repeated indices are summed up, unless specified otherwise. 

2.1. Lagrangian description 

Replacing v by d T X in Eq. (jl-ip . we can make the action in the Lagrangian 
description. The conservation laws of mass and entropy are given by 

pJ = Pinit and s = sjnit , (2-3) 

where p^t and sinit denote initial values of p and s, respectively. Thus, the action 
can be defined as 

rtfin r 
S L [p,s,X,K,X]= / dr d 3 a {JC(p,s,d T X)+K(pJ - p init ) + pAJ (s - s in it)} , 
Jt init Jv 

where K and A are undetermined multipliers introduced to keep the constraints 
(Eq. (|2-3[) ), In this variational calculus, both ends of a pathline are fixed, i.e., 

«Q(a,£ init ) = SXi{a,tsn) = , (2-5) 

where 5 indicates an infinitesimal variation. Let us write n for the unit normal vector 
of dV directed outside. The slip boundary condition, riid T Xi(a,T) = if a £ dV, 
means that a fluid particle that is initially in contact with dV remains in contact 
with dV although it can slip along the boundary. Thus, we have 

ni5Xi(a,T) = if a £ dV . (2-6) 

The stationary condition of Eq. (J2-4J) with respect to K, A, p, s, and Xi is 
respectively given by the two equations of Eq. (|2-3p , 

K = ~(d T Xif + h, (2-7) 

A = T , (2-8) 

where h is enthalpy defined as h = e+p/p. The undermined multipliers K and A are 
related to the physical quantities by means of Eqs. (|2-7p and (|2-8p . Note that surface 
integral terms, appearing when we applied the integration by parts in calculating 
Eq. (|2-9|) . vanish because of the boundary conditions (Eq. (|2-6|) ). 

From Eqs. (|2-7p - (|2-9p . we successfully obtain the Euler equation in the La- 
grangian description, 

d 2 x t _ dp dJ 

9 8t 2 d aj d(dXi/d aj ) ' l J 

The above is equivalent to the Euler equation in the Eulerian description, which is 
given by Eq. (|2-2ip below. We can show the equivalence by multiplying Eq. (|2-10|) 



and 



by J~ 1 (x,t), and replacing d 2 Xi/dr 2 and J~ 1 (dJ/d(dXi/da,j)) by (dt + v- V)i>j and 
dcij/dxi, respectively!^ As shown in Eq. (|2-5p . we fix both ends of a pathline. The 
Euler equation Eq. (|2-10|) is a second-order differential equation with respect to time 
r, and its solution has two constants of integration, which are determined by the 
fixed ends. 

2.2. Eulerian description 

In the Eulerian description, p, s, and v are regarded as functions of x = 
{xi,X2-,xs) and t, and we need not replace vi by d T Xi in Eq. (jl-ip . The bound- 
ary condition Eq. (|2-6|) gives 

n^(a;,£) = if x £ dV , (2-11) 

while the conservation laws of mass and entropy Eq. (j2-3|) can be rewritten into 

D t (p * 1) = , i.e., dtp = -V • (pw) , (2-12) 

D t s = , i.e., d t s = — v ■ Vs , (2-13) 

where Dt and *1 denote the Lagrangian time derivative equivalent to d T and volume 
element equivalent to J, respectively. See Appendix [A] for the details. If we follow 
a straightforward way from Eq. (|2-4)) . the action in the Eulerian description should 
be given by 

S E [p,s,v,K,X}= dt d 3 x lc(p, s,v)-J^+V ■ (pvyj-Xp(-^+v -Vsjl , 

(2-14) 
where k and A are undetermined multipliers. Keeping 

5p{x, t init ) = 8p(x, t fin ) = , (2-15) 

5s{x, t init ) = Ss(x, i fin ) = , (2-16) 

we find the stationary conditions of Eq. (|2-14p with respect to k, A, v, p, and s to 
be given respectively by Eqs. (|2-12p and (|2-13p . 

u = -Vk + AVs, (2-17) 

D t K = -^v 2 + h, (2-18) 

and 

D t X = T . (2-19) 

Note that surface integral terms appearing in the calculation vanish because of 
Eqs. (|2Tip . (|2-15p . and (|2T6p . As discussed later in ^4.2^ we have the same station- 
ary conditions, Eqs. (l2T2"1) . (l2 7 T3]l . and ^FWf - STWb . even if we replace Eqs. ([2T51) 
and (|2T6p by other restrictions of p and s at the initial and final times. Comparing 
Eqs. flMD and <^E> with Eqs (|2"T5j) and <^W), respectively, we find K = D t n and 
A = DtX. With the aid of Eq. (J2-13J1 . the Lagrangian time derivative of Eq. (|2-17p 
yields 

D t v = -VA« + (AA)Vs . (2-20) 



Substituting Eqs. (|2-18p and f)2-19[) into Eq. (|2-20|) . we obtain the Euler equation 

d 1 V?) 

—v + -Vv 2 -vx(Vxv) = . (2-21) 

at 2 p 

From Eq. (|2-17p . the vorticity is found to be given by 

u)eVxd = VAxVs. (2-22) 

Although the derivation of the Euler equation appears successful, the resultant vor- 
ticity Eq. (|2-22|) vanishes on the isentropic condition Vs = 0. It is the flaw mentioned 
in the introduction. 

§3. Fixing ends of a pathline 

Both ends of a pathline are fixed in the Lagrangian variational calculus, as shown 
in Eq. (|2-5p . It is thus natural that the ends are fixed in the Eulerian variational 
calculus, but not in ^2.21 

Let Ai (i = 1,2,3) be three scalar fields so that a pathline coincides with an 
intersection of hypersurfaces given by 

Ai(x,t) = constant for % = 1,2, 3 . (3'1) 

The determinant of the Jacobian matrix 

d(A 1 ,A 2 ,A 3 ) 



d(x 1 ,x 2 ,x 3 ) 



(VAi x VA 2 ) ■ VA 3 (3-2) 



gives a reciprocal of the volume element in the coordinates defined in terms of A. 
The mass conservation, the first equation of Eq. (|2-3p . is also represented by 

p (a(x, t), timt) (VAi x VA 2 ) • VA 3 = p{x, t) , (3-3) 

where A has variables x and t. 

We can take A for the Lagrangian coordinate a. Otherwise, we can consider A 
as an invertible function of a. The indefiniteness was also pointed out in Ref. [T0|) . 
Suppose, for example, we have 

A 1 = ai , A 2 = a 2 , and ^3 = 03 + /(01, o 2 ) , (3-4) 

where f(a\, a 2 ) is a two- variable function of a\ and a 2 . Because we have 

df df 

VA 3 = Va 3 + TT^Vai + ^Vaa , (3-5) 

oa\ oa 2 

we can define VA 3 so that it is normal to V-Ai(= V01) and Vj4i(= Va 2 ) by tuning 
the function /. If we take A to be an arbitrary invertible function of a, and substitute 
A, instead of a, into Eqs. (|3-4p and (|3-5p . we can also make Vvl3 normal to V.Ai(= 
V^4i) and V ' A 2 {= V^)- By using it, we can rewrite Eq. (|3-3|) into 

VA = P(x,t) ,„„, 

3 p(a(x,t),t init )(VA 1 xVA 2 )> { ' } 
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where A{ has variables x and t. 
Let us impose 

5A a (x,t init ) = 5A a (x,t &n ) = 0, for a = 1,2 . (3-7) 

Using Eqs. flEE]), ([321), and flj£I|), we have 

Ma(sMiiiit) = SA 3 (x,t &n ) = . (3-8) 

Thus, we can fix the ends of a pathline in the variational calculus by imposing 
Eqs. (J2-15P and f|3-T|) . Note that the values of Ai and A2 can be determined inde- 
pendent of the mass density p. 

Since v is tangent to the hyper surf aces, A\ and A2 must satisfy 

£Ma = 0,i.e., ^A a = -vVA a . (3-9) 

Hence, in the Eulerian description, we should minimize the action 

f tfin f V~-v 

Si[p,s,v,A,k,X,(3] = S-e[p,s,v,k,X]- dt dx 3 2Zp(3 Q D t A a , (3-10) 

Anit JV £= x 

with Eqs. (j2-15fl . (j2-16fl . and (|3-7p kept. Here, j3 a is the undetermined multiplier, and 
we write A and (3 for (A\,A2) and (/3i,/?2) respectively. The stationary conditions 
of Eq. (|3-10p with respect to A a are given by 

D t p a = , i.e., J^/3 Q = -v • V/3 Q . (3-11) 

Here, surface integration terms appearing in the integration by parts vanish because 
of Eqs. (|2-lip and (|3-7p . The stationary conditions with respect to v and p are 
respectively given by 



-V« + AVs + Y^ Pc^A a , (3-12) 



and 



a=l 

2 



D t K = ~v 2 + h + ^2f3 a D t A a , (3-13) 

which equals Eq. (|2-18|) because of Eq. (j39|) . The other stationary conditions with 

respect to k, A, p, and s are respectively given by Eqs. (|2T2p . (|2-13p . (|2T8p . and 

(|2-19|) . We obtain the Euler equation Eq. (|2-2ip from Eq. (|3T2p in the same way as 

we used in the preceding section. Values of 0% and fa are determined by the fixed 

ends, as discussed later in §4.21 From Eq. (|3-12p . the vorticity is found to be given 

by 

2 

u = VA x Vs + J2V/3 a x VA a . (3-14) 

a=l 

This term Yua=i ^ fia x V ' A a makes the flow rotational even on the isentropic con- 
dition Vs = 0, and the flaw mentioned in the introduction is removed. 



§4. Hamiltonian formulation 

Hamiltonian formulations in analytical mechanics can be regarded as a special 
case of more generalized formulation in control theory, known as Pontryagin's mini- 
mum principle .112) J2Q) T n £ |4,i[ we f] rs t give a brief review of this theory. In §4.21 we 
apply it to the dynamics of a perfect fluid to derive one of the canonical Hamilto- 
nian formulations with Clebsch potentials, which was previously derived in a different 

waj \mm 

4.1. Brief review of the control theory 

Let q represent the state of a system to be controlled, and u represent the input 
to this system, and we assume the time evolution of the state to be given in terms 
of a function of the state and input as 

j t q = F(q,u). (4-1) 

The optimal input is determined so that a cost functional 



f hn dtL(q(t),u(t)), (4-2) 

J iinit 



where L denotes a function, is minimized on condition that the initial and final states 
are fixed, i.e., 

<y«(*btt) = , (4-3) 

and 

5q(t &n ) = . (44) 

We define the undetermined multiplier, p, which is also called costate. The optimal 
input is obtained by minimizing 



S[q,P,u]= dt ll(q,u)+p- (—q-F(q, 



u 



m dt l-H(q,p,u)+p--q\ , (4-5) 

where H(q,p,u) is defined as 

H(q,p,u) = -L(q,u)+p-F(q,u) . (4-6) 

Let u*(q,p) denote the input minimizing Eq. (|4-5|) on condition that q and p 
are given, and it is necessary for u* to satisfy 

du* 

Introducing 

H*(q,p)=H(q,p,u*(q,p)), (4-8) 
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we define the preoptimized action as 

S*[q,p}= [ thn dt \-H*(q,p)+p.j t q\ , (4-9) 

which is not larger than S[q,p, u]. We can obtain the optimal input by solving 
the stationary conditions of Eq. (|4-9p with respect to p and q, which are given 
respectively by 

dt d Pi { ' 

and 

d Pi dH*(q,p) 



dt dqi 



(4-11) 



The boundary conditions for these equations are given by the initial and final states, 
which are assumed to be fixed in Eqs. (|4-3p and (|4-4p . If we impose not Eq. (|4-4p 
but a boundary condition 

p(tsn) = , (4-12) 

on the stationary conditions of Eq. (|4-9p . we can derive Eqs. (|4-10|) and (|4Tip . As 
discussed in Ref. |20|) . we can have more general conditions for the initial and final 
states. These conditions are called the transversality conditions in control theory. 
Let us assume q to represent the position of a material particle, and define F as 

F(q,u)=u, (4-13) 

and we find the formulation above to be equivalent to that of analytical mechanics 
for the particle. 

4.2. An application for a perfect fluid 

We can generalize the formulation in the preceding subsection to cases where 
the state variable and input variable are functions of space and time, respectively. 
Let us consider p, s, and A in ^3] to represent the state, and v to represent the input 
variable. We can identify Eq. flUO with a set of Eqs. (l2"T2l . (12T31) . and ([23]) by 
defining q = (p,s,A) and u = v. The costate is defined as p = (—K,—pA,—pP). 
The variations at the initial and final times satisfy Eqs. (|2-15p . (|216j) . and (|3-7p . 
The function corresponding to Eq. (|4-6p is given by 

2 

U(q,p,u) = -C + kV- (pv)+p\v-X7s + ^2p(3 a X7A a . (4-14) 



a=l 



Using Eqs. flH]) and flM!}, we find 



i'd 3 xH(q,p,u)= ld 3 x ln*(q,p)-^(v + VK-XVs-^2(3 a VA a ) 2 \ , (4-15) 



where we define 



U*(q,p) = p \ e(p, s) + ~(-Vk + XVs + ^ f3 a VA a ) 2 1 . (4-16) 

Thus, the velocity field given by Eq. ()3-12p turns out to minimize the action on 
condition that the state and costate are given. We find that Eq. (|3-10p satisfies 



f fln f { d 

Si[q,p,u]= / dt d 3 x \-H(q,p,u)+p- - 

•/tinit JV { « 

>j^dtjd 3 x L n *(q,p)+p-^\ , (4-17) 

which is the preoptimized action, S*[q,p\. Here, we have 



2 



dq dp ds ^ dA a 



a=l 



The stationary condition of the preoptimized action gives a set of Eqs. (|2-12jt . (|2-13|) . 
and (|3-9p . and a set of Eqs. (|2T8|) . (|2T9p . and (|3-lip . These sets can be rewritten 
respectively into 

dqi(x,t) _ dH*(q,p) y d dH*(q,p) 
dt dpi 4^i dxj d(dpi/dxj) 

and 

dpj(x,t) = &H*(q,p) | y d dH*(q,p) 
dt dqi 4^ dxj d(dqi/dxj) 

for i = 1,2,3,4, which correspond to Eqs. (|4-10p and (|4-llj) . respectively. We can 
rewrite Eqs. (|4-19j) and (|4-20p respectively into 



:! 



dH*(q,p) ,A d dH*{q,p) 



•(-,« = «(.,«+ I t jt \-^ l -'£^ w ^} (**> 



41 , ] &H*{q,p) A d &H*(q,p) 



f 1 , dH*{q,p) ^ 
Pi(x,ti)=pi{x,t- mit )+ / dt <{ h^ 



J dxj d(dqi/dxj) 



> .(4-22) 



In §0 we fixed the initial and final states by imposing Eqs. (|2-15p . (|2-16p . and (J3-7J) . 
These restrictions can be rewritten into Eqs. (|4-3|) and (|4-4|) . The initial value of the 
costate p(aj,ti n it), i-e., k(x , ti n n) , \(x , ti n n) , and (3{x,t- m \ t ), should be determined so 
that the given values of q(x,t m \t) and q(x,t^ n ) are compatible with Eqs. (|4-2ip and 
(14^21 . 
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We can relax the restrictions in the variational calculus. First, suppose that we 
do not impose Eq. (|3-7p for the final time. As discussed in relation to Eq. (|4T2p , we 
have /3(x,tft n ) = 0. By substituting it into Eq. (|3-11|) . we obtain 

(3(x,t) = 0. (4-23) 

Thus, the velocity field Eq. (|3-12p becomes irrotational on the isentropic condi- 
tion, like the velocity given by Eq. (|2T7p . Next, suppose we do not impose any of 
Eqs. ([EESD, d2HSD, and d22D for the final time, and we have Eq. (j¥2H and 

«(*,*&,) = A(a5, ton) = 0. (4-24) 

If we can solve Eqs. (|4-2ip and ()4-22p with their boundary conditions above, the 
solution uniquely determines p, s, A,k,X, and (3, and thus v, for any x and t con- 
sidered. 

§5. Discussion 

Various actions for a perfect fluid in the Eulerian description were proposed 
in previous works so that a rotational isentropic flow can be derived. Bateman 2 --' 
proposed an action with additional scalar field £ and r/ as 

S b [p, s, v, k, A, 77, C] = S E [p, s, v, k,X\- / dtdx 3 pQD t r] , (5-1) 

and derived 

v = - Vk + AVs + CVr? . (5-2) 

On the isentropic condition, the velocity field is given by 

v = -VK + CVr/. (5-3) 

The additional field £ and r/ are sometimes called Clebsch potentials because, apart 
from the variational principle, Clebsch 3 '' 4 ' had shown that a velocity field of a perfect 
fluid is given by Eq. (|5-3p . provided that the vorticity can be put in the form 

u> = VC x Vr] . (5-4) 

LirP 1 introduced more scalar fields to identify a fluid particle, mentioned the 
conservation law of identity without sufficient explanation of physical meanings of 
the identities and their conservation, and derived 

3 

v = -Vk + AVs + ^ftVAj , (5-5) 

j=l 

which can be rotational on the isentropic condition. He considered Ai to be a func- 
tion of the original position of the fluid particle. Later, Selinger and WhithanF* 
mentioned, "Lin's device still remains somewhat mysterious from a strictly mathe- 
matical viewpoint, but necessary for it seems to be firmly established as we proceed," 
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and also commented, "But it was noted by Clebsch that any flow velocity can be 
represented by" Eq. (|5-3p . This statement would be misleading because Eq. (|5-3p 
is valid if and only if Eq. (|5-4p holds™ Schutz-^ 1 generalized Bateman's variational 
principle for a relativistic perfect fluid, and mentioned, "The existence of an Eulerian 
variational principle may be a beginning. What is needed, I believe, is a Hamiltonian 
principle in a minimum number of variables. The present action principle seems to 
have 'too many' free variables". Yoshida ■ pointed out that an arbitrary velocity 
field can be given by Eq. (|3-12|) . not Eq. (|5-3p . in a way different from the way we 
used in $3j There, we first introduced Ai (i = 1,2,3), which can be the Lagrangian 
coordinates, but can be otherwise. Using the mass conservation Eq. (|3-6p . we next 
showed that we can fix both ends of a pathline without using A3. 

Clebsch potentials have been considered to be related with the Lagrangian co- 
ordinates JDJSJiEljIinHIIJ UMUB Selinger and WhitmarP mentioned that Clebsch po- 
tentials satisfy Eq. (|3-9p . Kambe^ considered Eq. (|55j) to give the transformation 
between the Lagrangian and Eulerian coordinates, and mentioned the indefinite- 
ness of Clebsch potentials. Yoshida 11 ' claimed that Clebsch potentials provide a 
one-to-one map between the Lagrangian densities in the Lagrangian and Eulerian 
descriptions. 

Although Clebsch potentials have been considered to label a fluid particle, why 
the label is required in the Eulerian variational calculus has not been clear in previous 
studies. We can fix both ends of a path line in the Eulerian variational calculus by 
imposing Eq. (|3-7p on the Clebsch potentials. In doing so, we can keep the same 
conditions in the variational calculus, as in the Lagrangian description. Thus, we 
can start with the same Lagrangian density to obtain the same velocity field, which 
can be irrotational on the isentropic condition, as in the Lagrangian description. As 
discussed in $H2J if we do not impose Eq. (j3-7|) . i.e., do not fix both ends of a path 
line, the resultant Eulerian variational calculus has no corresponding Lagrangian 
variational calculus, in which both ends of a path line are always fixed by Eq. (12 -5p . 

As referred to in Ref. [10]) . Lin's variational principle yields 

D t f3i = 0, for i = 1,2,3 (5-6) 

which means that /% is a function of the Lagrangian coordinates, i.e., Clebsch po- 
tentials Ai (i = 1, 2, 3), and the vorticity on the isentropic condition is given by 



3 



d(3i(A) 



w(Ai, A 2 , A 3 ) = Yl J ^ A - lvA ^ x VA * ■ ( 57 ) 

We showed that A3 is a function of p,A±, and A2 in Eq. (|3-6p . Thus, the values of 
/?1 and /?2 are given by functions of p, A\ , and A 2 because of Eq. (|3-11|) . Hence, the 
vorticity on the isentropic condition can be rewritten into 

u>(A 1 ,A 2 , P )=(^-^)vA 1 xVA 2 + Y,^VpxVA a . (5-8) 

Hamiltonian formulations for a perfect fluid also have been studied by some 
authorsP-OT^ 1 Arnold and KhesirASI studied a Hamiltonian formulation of an in- 
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compressible isentropic flow, and Morrison and Greene^* discovered an extremely 
complicated Eulerian Poisson bracket for an isentropic perfect fluid. Holm and Ku- 
pershmidt—' modified the Hamiltonian formulations^'''^'' by means of Clebsch po- 
tentials. As pointed out in Ref. llip . Casimir invariants, such as helicity, constrain the 
dynamics of a perfect fluid in these noncanonical Hamiltonian formulations Q2Hn) 
Note that the helicity is conserved only for a barotropic perfect fluid and an isen- 
tropic perfect fluid. On the other hand, as discussed in Appendix [Bj there is no 
Casimir invariant in a canonical Hamiltonian formulation. We obtain the same 
canonical Hamiltonian formulation as derived in Ref. 112]) by regarding the velocity 
as the input variable. 

We have assumed fluid dynamics in a fixed container. This assumption can be 
relaxed; our discussion remains valid when we consider the dynamics of a perfect 
fluid in an infinitely large spacial region and assume that the mass density tends to 
become zero far from a certain point. Our discussion can be applied to a relativistic 
prefect fluid, where we should start with the Lagrangian density slightly different 
from Eq. (jl-ip . We show the details in Appendix O for completeness of our study. 
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Appendix A 

Lagrangian Time Derivative 



The Lagrangian time derivative d T = Dt is given by dt + v ■ V for a scalar, 
dt + V(v- ) — v x Vx for a cotangent vector, dt — V x (»x ) + t>(V- ) for an axial 
vector, and dt + V • (v ) for the volume. These expressions are unified as 



D t = d t + L v , (A-l) 



where L v denotes the Lie derivative along the vector field v. Note that Dt is commu- 
tative with exterior derivative d, and thus gradient V- and rotation Vx . Because 
p and s are 0- forms, and the volume element *1 is a 3-formpQ the conservation laws 
of mass and entropy are respectively found to be given by 



D t (p * 1) = and D t s = . (A-2) 



13 

Appendix B 

Canonical Poisson Bracket 

Let / and g denote functions of q and p, which are TV-dimensional vectors, and 
we define their Poisson bracket as 



{/,<?} = £ 



5qi Spi 5pi 5qi 



where we define 



i=l 



5 d v^ d d 



(B-l) 



— - = — -- V— - — — -— — (B-2) 

<% dqi ^ dxj d{dqi/dxj) 



3 

5 d A 



fyi dpi ^ i dxjd(dp i /dx j ) 



We can rewrite Eqs. ()4-10p and ()4-lip into 

j t q = {q,H*(q,p)}, (B-4) 

j t P = {p,H*(q,p)} , (B-5) 

respectively. This Hamiltonian formulation is canonical because the associated Pois- 
son bracket can be rewritten by means of the symplectic matrix, i.e., 

ySqi' 5piJ \ -I J ySqi' SpiJ 

where I and the superscript t are the N x N unit matrix and the transposition, 
respectively. A Casimir invariant C is defined so that {F, C} = for any F, and is 
a conserved quantity since dC/dt = {H, C} = 0. Thus, if C is a Casimir invariant, 
we have 

I \ (5C 5C X ' 



i oaS'KJ =0> (B7) 

which tells that C does not depend on any of q and p, i.e., that C is trivial. 

Because Eqs. (|4T9p and (|4-20p can be rewritten by means of the symplectic 
matrix, our Hamiltonian formulation given in ^4.21 is canonical. 

Appendix C 

Formulation for a Relativistic Perfect Fluid 



Let us consider a relativistic perfect fluid in a four-dimensional space-time region 
Q. We redefine p and s as the particle number density and the entropy per particle 
in the local rest frame of matter, respectively, and define u as the four-velocity field. 
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We assume that metric tensor g( , ) is given by other materials, and put the speed 
of light equal to unity. The normalization of the four- velocity field u is given by 

#(u, it) + 1 = 0. (C-l) 

The constraints Eqs. (|2-12p . (|2-13p . and (|3-9p are respectively rewritten as 

L u (p * 1) = 0, L u s = 0, and L u A a = 0, for a = 1, 2 . (C-2) 

The Lagrangian density for a relativistic perfect fluid is given by 

C(p,s) = -pe. (C-3) 

Let *1 denote the four-dimensional volume element, and we can modify the action 
Eq. (|5TD|I as 

S(p,s,A,K,\,/3,u,~/)= |*LC0,s)+7* l{g(u,u) + 1} 

2 

-kL u (p* l)-p\* lL u s-^2p(3 a * lL u A a j , (C-4) 

where 7 is an undetermined multiplier for the normalization. We still have Eq. (|2-lip . 
and impose restrictions Eqs. (|2-15|) . (|2-16|) . and (|3-7|) in the variational calculus. The 
stationary conditions corresponding to Eqs. (|2-18p . and (|2-19[) are respectively given 
by 

-h + L u k = , (C-5) 

-T + L u X = 0, (C-6) 

L u Pa = , (C-7) 

2 
2-/g(u, )+dn- Xds-^2 P^dA a = , (C-8) 

a=l 

where enthalpy h, pressure p, and temperature T are defined in the local rest frame 
of matter. From Eqs. flgEQ) , (RT2]) . (RT31> - (l(T8l> . we obtain 27 = h and thus have the 
equation for the four-momentum field 

2 
hg(u, ) + dn- Xds - V] /3 a d^ a = , (C-9) 

i=l 

of which the Lie derivative gives the Euler equation 

L u {hg(u, )} + dp/p = . (C-10) 

Defining 

v= — ,— ,— , (C-ll) 

V^O «o MO/ 
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we have no = l/\/l — V 2 from Eq. (jC-ip . As in §4.2t we take p = puo, s, and A to 
be the state q, take — k, —p\ and — p(3 to be the costate p, and take v to be the 
input. Because Eq. (JC-2P gives 



^ = -v • 0*0 

5t s 



-w • Vs , 
-v ■ X7A a , 
we can find the function corresponding to Eq. (|4- 14j) to be given by 



dt 



(C-12) 
(C-13) 

(C-14) 



H{p, s, A, A, k, (3, v) = -£(p\/l - v 2 , s) + kV ■ (pv) + pXv ■ Vs + ^ p/3 a i> • VA a , 

(C-15) 
and the optimized input v*(q,p) to satisfy Eq. (|C-9p . 
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